Assortative Mixing by Degree Makes a Network More Unstable 
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We investigate the role of degree correlation among nodes on the stability of complex networks, 
by studying spectral properties of randomly weighted matrices constructed from directed Erdos- 
Renyi and scale-free random graph models. We focus on the behaviour of the largest real part of 
the eigenvalues, A max , that governs the growth rate of perturbations about an equilibrium (and 
hence, determines stability). We find that assortative mixing by degree, where nodes with many 
links connect preferentially to other nodes with many links, reduces the stability of networks. In 
particular, for sparse scale-free networks with N nodes, A max scales as N a for highly assortative 
networks, while for disassortative graphs, A max scales logarithmically with N. This difference may 
be a possible reason for the prevalence of disassortative networks in nature. 

PACS numbers: 89.75.Hc, 05.65.+b, 05.45.-a 



Whether complex systems are more stable than sim- 
pler ones has been a contentious issue for a long time 
pJQ- Stability of a N~ dimensional dynamical system 
x = —x + F(x), (e.g., a network of N nodes, each node 
i associated with a variable Xi whose time evolution is 
governed by a general nonlinear function F) is measured 
by the rate at which a small perturbation about an equi- 
librium state x* decays with time. This is determined 
by the largest real part, A max , of the eigenvalues of the 



Jacobian matrix Ji, 



dFi /dxj 



representing 

the interactions within the system. If A max is positive, 
the effect of a small perturbation will grow exponentially 
with time, and the system will be rapidly dislodged from 
the equilibrium state. In the absence of detailed knowl- 
edge about the interactions among the constituent ele- 
ments of a complex system which determine Jij , ensem- 
bles of randomly constructed matrices have been stud- 
ied 3]. Results from several studies of random matri- 
ces, symmetric as well as asymmetric |J, |5(, show that 
Amax ~ V NCa 2 where N is the system size, C is the 
connectivity (the probability that there is a link between 
any pair of nodes in the network) and a is a measure 
for interaction strength between the nodes. Assuming 
self regulation, e.g., J a — — 1, it follows that a transition 
from stability to instability takes place as the network ei- 
ther grows in size or becomes more densely (or strongly) 
connected. This result, implying complexity promotes 
instability, has been shown to be remarkably robust with 
respect to various generalisations [E 0> HI • 

There is a direct correspondence between the nature 
of the matrices J and the structure of the underlying di- 
rected network, because two nodes i,j are connected iff 
Jij =/= 0. Clearly, the random matrices investigated in 
these previous studies correspond to Erdos-Renyi (ER) 
random graphs 0. Recently, it has been pointed out in 
a series of studies that many real world networks differ 
from these purely random graphs 0,0, e.g., in having 
highly skewed degree distributions that typically follow a 



power law with the exponent 7 between —2 and —3 |l2| . 
Previous studies have investigated spectral properties of 
undirected binary graphs with power law degree distri- 
bution 0,0,^3- However, there are further detailed 
features that make real world networks appear funda- 
mentally different from random graphs. In particular, 
in almost all biological networks, nodes of high degree 
tend to avoid being connected to other highly connected 
nodes, i.e. these networks show disassortative mixing 
|l6l Il7| . Understanding why networks exhibit these pat- 
terns is one of the challenges of current research ^| . 

In this paper we show that assortativity can signif- 
icantly influence network stability. For this, we study 
spectral properties of randomly weighted matrix ensem- 
bles into which we introduce correlations between the 
matrix elements Jij that reflect graph properties found 
in real world networks. We first examine the growth of 
Amax with network size in matrix ensembles correspond- 
ing to graphs with Poisson degree distributions. Increas- 
ing the assortativity of such graphs results in increas- 
ing the largest eigenvalue, implying decrease in stability. 
Next, we study the behaviour of growing networks with 
scale free degree distribution. We observe the eigenvalue 
distribution (that departs significantly from the spectra 
of the corresponding adjacency matrices) and determine 
the scaling of A ma x with network size. Again we find that 
assortative mixing leads to a loss of stability. Our results 
seem to imply that the disassortative mixing seen in bio- 
logical networks could be a consequence of evolutionary 
optimisation of these networks to minimise the effect of 
dynamical fluctuations. 



Every directed graph can be represented by its adja- 
cency matrix Ajj , where Aij = 1 if there is a link from 
j to i; = otherwise. Following May 0, we generate 
the matrix Jij from Aij of a network with desired prop- 
erties, by replacing the non-zero elements of the latter 
matrix with entries drawn from a Gaussian distribution 
having mean zero and variance a 2 . Note that, as the re- 



2 



-0.12 
-0.14 
-0.16 



<< -0.18 



-0.2 
-0.22 




0.5 1 1.5 2 GD-, 



0.05 



-0.6 -0.4 -0.2 
assortativeness a 



0.2 



FIG. 1: Dependence of the average of largest real part of 
the eigenvalues, (A max ), on assortativity for Poisson degree 
distributed networks of size ./V = 100. Data points represent 
averages over 1000 independent runs. The inset shows the 
tail of the real axis projection of the spectrum for random 
matrices derived from non-assortative (a ~ 0, filled circles) 
and strongly assortative (a — 0.8, open circles) networks. The 
tail becomes more extended for strongly assortative networks. 

sultant matrices are not symmetric, the eigenvalues are 
expected to be found in the complex plane. In this ap- 
proach, correlations between the matrix elements are de- 
termined by the structure of the underlying graphs. In 
particular, third and higher order correlations are related 
to the network assortativity, 

r(EtiJ^)-(lEta0i + ^)) 2 
r(Eti -M + *?) - (i(Ef=i hV< + k *W ' 

where ji,ki are the degrees of nodes connected by the 
z-th link and L ~ N(N — 1)C is the total number of links 
in the network. 

First, we consider matrix ensembles constructed from 
graphs with Poisson degree distributions. We start with 
an ER graph (for which a w 0) and then introduce as- 
sortativity by applying the method proposed in Ref. . 
The algorithm essentially consists of performing a series 
of rewiring procedures, such that, to obtain an assorta- 
tive network, the steps leading to higher (lower) assor- 
tativity are accepted with increased (decreased) proba- 
bility. Fig. n illustrates the dependence of (A max ) on 
the assortativity. These data show a clear trend towards 
larger (A max ), i.e. less stability, as networks become more 
assortative. Note that, changing the assortativity by this 
rewiring procedure also increases the clustering coeffi- 
cient, c. To correct for this, we use the same scheme to 
change the clustering coefficient, while keeping the degree 
distribution fixed. As this leaves the assortativity un- 
changed, we can generate networks with a specific value 
of c while keeping a constant. In this case, there is no 
statistically significant change of (A max ) with c (data not 
shown) . Combining this result with Fig. ^ allows us to 
conclude that high assortativity leads to lower stability 




FIG. 2: (Top) Real axis projection of the spectral densities 
(shifted to the right by AA = 1) for matrices constructed 
from BA networks (filled circles) and randomised BA net- 
works (open circles) with N = 5000, m — 10, and a — 0.2. 
The left inset illustrates the exponential tail of the spectrum, 
while the right inset shows that the spectrum follows a power 
law -P(A) ~ A~ 1//2 for small A in the BA case, but not after 
randomisation. The comparison illustrates that the 'dynamic 
correlations' in the BA network have an effect on the spec- 
trum. The data represent averages over 100 independent runs. 

for networks with Poisson degree distribution. 

Next, we look at random matrices constructed from 
scale-free graph substrates. The networks are generated 
by adding at each time step a new node with m links 
that connect to existing nodes following the Barabasi- 
Albert (BA) preferential scheme presented in Ref. [l^ . 
To avoid 'dynamic correlations' resulting from the BA 
procedure, we subsequently randomise the networks us- 
ing the link swapping technique discussed in Ref. |2fj| . 
Corresponding Jij matrices are constructed as before. In 
Fig. |21 a comparison of numerical data for the spectral 
densities of unrandomized and randomised BA networks 
demonstrates that the dynamical correlations in the for- 
mer indeed have an impact on the spectral density. While 
both spectral densities have exponentially decaying tails 
[Fig. |21 (left inset)], only the spectral density for the ma- 
trix ensemble constructed from unrandomized BA net- 
works seems to obey a power law P(X) ~ A -1 / 2 for small 
A [Fig. [21 (right inset)]. Note that, this result obtained for 
weighted, directed graphs differs significantly from the 
eigenvalue spectrum of adjacency matrices of scale-free 
networks, which exhibit power-law tails and exponential 
decay (or "triangle- like" shape) at the centre [l3|, |lj] . 

We now observe the scaling of (A max ) with network 
size N. The well known results for ER graphs suggest 
that A max essentially is a function of the number of links 
in the network and their average link strength. Hence, 
to explore the dependence of (A max ) on network size, we 
look at two cases. 

In the first case, as seems realistic for, e.g., food webs 
[2l| . we keep the connectivity C of the networks con- 
stant as we increase the network size N. It follows that 
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FIG. 3: Scaling of (Amax) with TV for exponential (open 
squares), BA (filled circles), and Bianconi (open circles) net- 
works with constant connectivity C = 0.02. The line indicates 
a power law with exponent 1/2, implying that, in the asymp- 
totic limit, A m ax ~ TV 1//2 for all three types of networks. 

the number of links per node added, m, becomes a lin- 
ear function of network size. Fig. [3] shows the scaling 
of the average of largest real part of the eigenvalues, 
(A m ax), with network size TV. For three types of grow- 
ing networks, viz., having exponential (i.e., no preferen- 
tial attachment) and scale-free degree distribution with 
exponents 7 = — 3 (the original BA preferential attach- 
ment algorithm) and 7 rs —2.1 (Bianconi networks, con- 
structed after a modification of the BA procedure follow- 
ing Ref. we find the same scaling behaviour as in 
ER graphs, A max ~ TV 1 / 2 , which contrasts with the rela- 
tion Ai ~ TV 1 / 4 reported for adjacency matrices of large 
undirected BA networks Variation among different 
network topologies are chiefly expressed in different finite 
size behaviour and proportionality constants. 

In general, it seems that networks with degree dis- 
tributions having a less extended tail are more stable. 
This appears to be supported by Gersgorins theorem |23| 



which gives the bound: A max < Ru 



1, where R n 



gree, k n 



is related essentially to the largest de- 
of the network. From TVy^2° fc -7 < 1, 



where k are node degrees (see e.g., Ref. |2J|), one obtains 
the following dependence of fc max on TV for the differ- 
ent networks: fc max ~ TV - 91 (for the Bianconi network), 
Lax ~ TV 1 / 2 (for the BA network) and fc max ~ log(TV) 
(for the exponential network). This ordering according 
to decreasing growth rates of fc max with TV agrees very 
well with the ordering of the network types according to 
decreasing stability in Fig. 

As already mentioned, in the preceding case we con- 
sidered the network connectivity C to be constant as TV 
increases. However, most real networks are sparse, i.e. 
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FIG. 4: (Top) Scaling of (A max ) with network size N when 
connectivity C decreases with increasing N. Data for net- 
works constructed from the BA procedure (filled circles) are 
compared with networks optimised for assortativity (open 
boxes, a m 0.24, c w 0.45) and for disassortativity (open 
diamonds, a ~ —0.2, c = 0). For the assortative networks, 
Amax ~ iV ' 2 (dashed line), while for disassortative and uncor- 
rected BA networks, A ma x ~ log(A r °' 12 ) (solid line) for large 
N. (Bottom) Scaling of the excess density / rca i of real eigen- 
values with TV: / rea i oc TV -0 ' 1 (dashed line) for the assortative 
networks, and / rea i oc N~ ' 45 (solid line) for uncorrelated BA 
and disassortative networks. Data points represent averages 
over 1000 runs (N < 800) and 200 runs (TV > 800). 

C - TV" 1 (see, e.g., Ref. H3). Therefore, in the sec- 
ond case, we measure (A max ) in networks with increasing 
size, while keeping m constant (resulting in decreasing 
connectivity C with increasing TV) . We now compare the 
behaviour of A max for assortative and disassortative net- 
works, constructed using the same method as before. For 
the results reported in this paper, we choose m = 2, for 
which we verify numerically that, after randomisation, 
the resulting networks are still connected. Optimising 
for a high assortativity, however, tends to fracture this 
component, e.g., we find that networks with a « 0.25 typ- 
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ically have a giant component that comprises only 80% of 
the links in the network. Note that, this fraction of links 
comprising the largest connected cluster, is independent 
of the number of nodes N. 

Fig. 21 (top) shows the scaling of (A max ) with N, for 
constant m. As in the case of the ER graph-based matri- 
ces, we again find that high assortativity leads to lower 
stability. For the assortative scale- free networks with a « 
0.25 we find that A max ~ TV 2 , whereas for the uncorre- 
lated and disassortative networks, A max ~ log(iV ) for 
large N. This difference in the scaling behaviour is sug- 
gestive of a huge difference in stability for large network 
sizes. 

The reason for the different scaling of A max in assor- 
tative and disassortative networks seems to be found in 
different scaling behaviours of the excess densities / rea i 
of real eigenvalues [Fig. 0] (bottom)] . For random ma- 
trices based on a ER graph topology, it is known that 
/ rca i oc N^ 1 / 2 for large N 0. For uncorrelated BA net- 
works we find that / rea i oc N~ s , 8 = 0.45 ±0.05, which 
is very close to the value for ER graphs. However, for 
assortative networks, we observe S = 0.1 ± 0.05, which 
results in a much more gradual decline than in disassorta- 
tive networks. Hence, for the same network size N, more 
eigenvalues lie on the real axis in assortative networks 
and thus the probability of having higher maximum real 
parts is increased. We also note that, as in the case of 
matrices based on graphs with Poisson degree distribu- 
tion [Fig. [T] (inset)], the tails at the spectrum edge show 
a more gradual decay for assortative networks, compared 
to uncorrelated and disassortative graphs (not shown). 

In this paper, we have presented an analysis of ran- 
domly weighted matrices whose corresponding networks 
are constructed by starting from well known graph mod- 
els such as ER and scale-free networks, and then using 
a link swapping algorithm to generate degree-correlated 
network ensembles with a desired assortativity. As the 
most important result of the paper, we find that assor- 
tative mixing by degree tends to severely destabilise net- 
works. While the largest real part of the eigenvalues 
scales as a power of the network size in highly assorta- 
tive networks, in the disassortative case we only observe 
a logarithmic dependence. 

The question has been raised why almost all biologi- 
cal networks are disassortatively mixed 0] • It nas been 
suggested that disassortativity has its origin in the re- 
striction that no two nodes have more than one link be- 
tween them. However, a recent study |25j shows that 
while this mechanism does indeed give rise to the kind of 
correlations actually observed in certain networks (e.g., 
the Internet), only a part of all the measured correla- 
tions can be accounted in this way. Other studies have 
pointed out that assortative mixing has a big impact on 
the percolation behaviour of networks 0]. This implies 



that disassortative networks might have been favoured 
because perturbations have a reduced chance to propa- 
gate. Our results add another viewpoint to this issue. 
The results reported here imply that disassortative net- 
works are more resistant to the effect of dynamical fluctu- 
ations than assortative networks j2(| . One may therefore 
speculate that, an evolutionary drive towards systems 
with reduced intensity of fluctuation (homeostasis) might 
be one of the reasons for the prevalence of disassortative 
networks in the biological world today. 
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